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Abstract. We present an asymptotic form of the Boussi- 
nesq equations in the limit of small Peclet numbers i.e. 
when the time scale of motions is much larger than the 
time scale of thermal diffusion. We find that, in this 
limit, the effects of thermal diffusion and stable strati- 
fication combine in a single physical process. This pro- 
cess is an anisotropic dissipation (not effective for horizon- 
tal motions) which acts primarily on large scale motions. 
The small-Peclet-number approximation presents also the 
great practical interest to avoid the numerical difficulty 
induced by the huge separation between the diffusive and 
dynamical time scales. The relevance of this approxima- 
tion to study the flow dynamics within the stellar radiative 
zones is considered. 
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1. Introduction 

The understanding of the flow dynamics within stellar ra- 
diative zones constitutes a major challenge for the cur- 
rent theory of stellar evolution. These motions transport 
chemical elements and it turns out that their contribu- 
tion might reconcile the existing models of stellar structure 
with t he obs ervations of the surface abundances (Pinson- 
neault 1998 ). By transporting angular momentum, such 



flows also play an important role in the evolution of star's 
rotation. In particular, they could explain the nearly solid 
body rotation of the solar radiative zone which has been 
revealed by helioseismology (Gough et al. 1996). 

We consider here the effects of the very high thermal 
diffusivity of stellar interiors on the dynamics of these mo- 
tions. In most cases, radiation dominates the thermal ex- 
changes within radiative zones. This heat transport is so 
efficient that the thermal diffusivities associated with the 
radiative flux are larger by several orders of magnitude 
than the thermal diffusivities encountered in colder me- 
dia like planetary atmospheres. For example, the thermal 
diffusivity varies between 10 5 and 10 7 cm 2 s _1 inside the 



sun whereas it is equal to 0.18 cm 2 s 1 in the standard 
conditions of the terrestrial atmosphere. 

This property of the stellar fluid is expected to strongly 
affect the flow dynamics especially inside the stably strat- 
ified radiative zone where the time scale of thermal diffu- 
sion appears to be shorter than the dynamical time scale 
characterizing radial motions. Helioseismology data show 
that the thermal structure of this region is very close to 
the one predicted by hydrostatic models, indicating that 
existing fluid motions are not fast enough to modify sig- 
nificantly the thermal structure built up by the radiative 



flux (Canuto & Christensen-Dalsgaard 1998). 

Qualitatively, the damping of temperature fluctuations 
by thermal diffusion is expected to have two main effects 
on the dynamics. The first one is to reduce the ampli- 
tude of the buoyancy force. This restoring force acts on 
fluid parcels displaced from their equilibrium level and is 
proportional to the density difference between the parcel 
and its environment. Since density fluctuations are pro- 
portional to temperature fluctuations for incompressible 
motions, fast thermal exchanges reduce the force ampli- 
tude. An important consequence of this effect is to favour 
the onset of shear layer instabilities in stably stratified 



layers (Dudis|l974| Zahn|l97J). 

The second main effect of the thermal diffusion is to 
increase the dissipation of kinetic energy. Any vertical mo- 
tion in a quiescent atmosphere induces a work of the buoy- 
ancy force so that a fraction of the injected kinetic en- 
ergy is necessarily transformed into potential energy. If 
the fluid parcels could "fall" adiabatically towards their 
equilibrium position, all the stored potential energy could 
return back to kinetic energy. However, the damping of 
temperature fluctuations provokes an irreversible loss of 
kinetic energy. A simple example of this process is the 
damping of gravity waves. 

Both effects of the thermal diffusivity are thus op- 
posed. While a decrease of the buoyancy force amplitude 
reduces the associated work and thus the amount of ki- 
netic energy extracted, the second effect increases the 
fraction of the kinetic energy which is irreversibly lost. 
Then, for a given mechanical forcing, a relevant question 
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is whether a larger thermal diffusivity reduces or enhances 
the kinetic energy of the flow. 

We are lacking quantitative results especially in non- 
linear regimes to answer such a basic question and more 
generally to understand the effect of thermal diffusivity in 
a stellar context. This situation is partly due to the dif- 
ficulty to reproduce flows with realistic Prandtl numbers 
either in laboratory experiments or by numerical simu- 
lations. The Prandtl number P r = v/k which compares 
the kinematic viscosity v and the thermal diffusivity k, 
varies between 10 -6 and 10~ 9 within the sun whereas it 
is equal to 0.7 in the air. Although some fluid like metal 
liquid may have small Prandtl numbers in laboratory con- 
ditions (P r = 0.025 for mercury, see for example Cioni et 
al. 1997)), these values remains far from the stellar case. 
The severe numerical limitation is explained by the huge 
separation between the time scales of viscous dissipation 
and thermal diffusion. The computation of both processes 
over a few dynamical times would require a prohibitive 
amount of computer time. 

In this paper, we investigate the limit where the time 
scale characterizing the thermal exchanges is much shorter 
than the time scale of the motions (the ratio between 
both time scales defines the Peclet number). In Sect. 2, 
an asymptotic form of the governing equations is derived 
in the context of the Boussinesq approximation. Evidences 
that these asymptotic equations actually approximate the 
Boussinesq equations for small Peclet numbers are pre- 
sented in Sect. 3. Then, in Sect. 4, the elementary prop- 
erties of the small-Peclet-number equations are described, 
emphasizing their theoretical and practical interests. Fi- 
nally, the relevance of this approximation in a stellar con- 
text is commented in Sect. 5. 



09 1 
— +u-V6 + w = — V 2 £ 

at P e 



V • u = 0, 



(2) 
(3) 



where, u = ue x + ve y + we z is the velocity vector, p the 
pressure and 9{x,y,z) = T{x,y,z) — T l {z) the tempera- 
ture deviation from the initial temperature profile. The 
z axis refers to the vertical direction, while the x and y 
axis refer to the horizontal ones. In the heat equation, the 
third term of the left hand side corresponds to the vertical 
advection of temperature against the mean temperature 
gradient dT l {z)/dz. This gradient is equal to unity in the 
dimensionless unit. To non-dimensionalize the equations 
we used the velocity scale [/*, the length scale L„, the dy- 
namical time scale t^ = L*/U*, the pressure scale QoU* 
and the temperature variation AT*. 

The system is then governed by the Richardson num- 
ber, the Peclet number, P e , and the Reynolds number, 
R e , respectively defined as 



P 



The Richardson number is the square of the ratio between 
the dynamical time scale £d and the buoyancy time scale 
<b = 1/-/V*. The thermal diffusivity n appears in the Peclet 
number which compares the thermal diffusion time scale, 
t K = L%/k with the dynamical time scale. The Reynolds 
number is the ratio between the viscous time scale L%/v 
and the dynamical time scale. 

In the limit of small Peclet number, we assume that 
the solutions u and 9 of the Boussinesq equations behave 
like Taylor series: 



2. Derivation of the small-Peclet-number 
approximation 

We restrict ourselves to a fluid layer embedded in an uni- 
form vertical gravity field and bounded by two horizontal 
plates. A mechanical forcing is assumed to drive motions 
which can be described by the Boussinesq approximation. 
We do not need to specify the forcing for the moment, we 
only assume that it introduces a velocity scale U*. The 
temperature is fixed on both plates so that a linear dif- 
fusive profile denoted T l (z) is established initially. The 
dynamical effect of the stable stratification is measured 

1/2 

by the Brunt- Vaisala frequency, — (fig AT*/ L») ' , 
where g denotes the gravity acceleration, [3 is the thermal 
expansion coefficient, AT* the temperature difference be- 
tween the upper and lower plates and L* the distance 
separating the plates. 

In the context of the Boussinesq approximation, the 
governing non-dimensional equations read: 



U = U + P e Ul + K U 2 



= 9 + P t 



(4) 
(5) 



Note that in the context of the Boussinesq equation, the 
pressure is an intermediate variable determined by the in- 
compressibility condition (j^). By inserting these asymp- 
totic expansions in the heat equation, we find at the zero 
order in P P : 



V 2 6» n = 0. 



(6) 



Since the temperature remains fixed to its initial value on 
both bounding plates, temperature deviations vanish on 
both plates. Then, Eq. (ph implies 



= 0. 



(7) 



Thus, at the lowest order in P e , the Boussinesq equations 
reduce to the Navier-Stokes equation: 



du 



1 9 

-Vp + R l 9e z + — V 2 u, 

tie 



(1) ^£ 

1 ' dt 



u • Vu 



1 9 



(8) 
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together with the incompressibility condition, 

V ■ u = 0. (9) 

At this order, the dynamical and thermal equations are 
decoupled. The coupling is recovered at the first order in 

-^+uo-Vui+ui-Vu = -Vpi+R i d 1 e z + — V 2 ui,(10) 

Ot H e 



w =V 2 8 1 , 
V ■ ui = 0. 



(11) 
(12) 



Solutions u = Uo + P e Ui, 9 = 6o + Pe@i valid up to the first 
order in P e must satisfy the above system of equations (Q) , 
(|, (|, ©, (0), (|ll). 

We note that the Lagrangian derivative of temperature 
deviations does not appear in the heat equation of this 
system. Thus, at the first order in P e , one would have 
found the same system of equations for Uo, Ui, Qq, 6\ 
if the Taylor series had been introduced in the following 
equations: 



du 



u • Vu = -Vp + R t 9e z 



Rp. 



P P w = V 



V-u = 0. 



(13) 

(14) 
(15) 



Therefore, if u and 9 actually behave as Taylor series for 
small Peclet numbers, the solution of the above equations 
is identical to the solution of the Boussinesq equations up 
to the first order in P e . 

The unique difference with the Boussinesq equations 
comes from the heat equation. Physically, the process lead- 
ing to the balance P e w = V 2 6* can be described as follows: 
For large values of the thermal diffusivity, the temperature 
fluctuations are expected to be small and the mean tem- 
perature stratification to remain unchanged by the me- 
chanical heat flux. However, vertical motions advecting 
fluid parcels against the mean temperature gradient al- 
ways produce temperature deviations and, unlike the non- 
linear advection term u.V#, this generation process does 
not depend on the amplitude of the temperature devia- 
tions. As fluid parcels go up (or down) in a mean tem- 
perature gradient, the amplitude of the temperature de- 
viations tends to increase continuously. In the mean time, 
thermal diffusion tends to reduce these temperature de- 
viations. Inspection of the heat Eq. (||) shows that this 
diffusive process can lead to a stationary solution, namely 
P e w — V 2 9. Clearly, if the time scale of the vertical mo- 
tions is very slow compared to the diffusive time scale, one 
expects that this stationary solution is practically instan- 
taneously reached. Again, it describes a balance between 
thermal diffusion and vertical advection against the mean 
temperature stratification. 



In the remainder of this paper, we will refer to the set 
of equations (pi), (|l4|), (15), as the small-Peclet-number 



equations or as the small-Peclet-number approximation. 
However, formal mathematical proof that u and 9 actually 
behave as Taylor series does not exist in the general case. 
Then, to prove that the small-Peclet-number equations ac- 
tually approximate the Boussinesq equations in the limit 
of small Peclet number, specific cases have to be consid- 
ered. In the next section, we shall present two types of 
linear flows where the validity of the small-Peclet-number 
approximation can be proved. Some evidences will also 
be given for a non-linear flow. The theoretical and practi- 
cal interest of the small-Peclet-number equations will be 
emphasized in Sect. 4. 

3. Validity of the small-Peclet-number 
appr oximat ion 

The first example we consider is that of small amplitude 
perturbations in a linearly stably stratified atmosphere. 
The perturbations are resolved into modes proportional 
to exp(— at)exp[i(k x x + k y y + k z z)\ where a is a complex 
number and k x , k y , k z , represent the horizontal and verti- 
cal wave numbers of the perturbation. In the following, the 
dispersion relation obtained using the Boussinesq equa- 
tions is compared to that derived from the small-Peclet- 
number equations. 

The calculation is conducted for two dimensional dis- 
turbances (k y = 0), but the three-dimensional case can be 
readily recovered replacing k 2 by k 2 + k 2 in the follow- 
ing expressions. To simplify the presentation we also limit 
ourselves to the inviscid case. It has been verified that our 
conclusions are not affected by taking into account the 
viscosity. 

Using the Boussinesq equations, the dispersion relation 

is: 



a — <jtv + Og = 

whereas the dispersion relation reduces to 



(j j 1 



(16) 



(17) 



in the context of the small-Peclet-number equations. In 
these expressions, 



k 7 



a 'I 



P 



is the damping rate associated with a pure thermal diffu- 
sion, and 



cb — v Ri 



k x 



is the frequency of gravity waves in absence of diffusive 
processes. 

We observe that, in the context of the small-Peclet- 
number approximation, all disturbances are damped with 
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a rate equal to o 2 b /ot- On the contrary, the dispersion re- 
lation of the Boussinesq equations shows different types of 
solutions. These solutions are now analyzed for increasing 
values of the thermal diffusivity. 

At small thermal diffusivity, solutions of the dispersion 
relation correspond to gravity waves damped by thermal 
diffusion. The two roots of Eq. ( [l6|) correspond to two 
gravity waves propagating in opposed direction. Increasing 
the thermal diffusion reduces the wave frequency until the 
roots of ([If]) become purely real and the associated modes 
damped without propagating. This occurs when 



Pe < 



3/2 



2 \/ jRij hi f 



It is important to note that, as long as k z is not equal 
to zero, there always exists a Peclet number such that 
this expression is verified for all k x and k z . If this was 
not the case, one could have gravity waves whatever the 
value of P e . Then, the small-Peclet-number approximation 
would not be valid since gravity waves are absent in this 
approximation. 

In deriving the small-Peclet-number equations, we re- 
stricted ourselves to a fluid layer bounded vertically. Ver- 
tical wave numbers have therefore a lower limit, fc™ m , 
so that all modes are damped without propagation if 

Then, the two distinct roots of the dispersion relation 
correspond to two damping modes. By further increasing 
the diffusivity, the damping rates take increasingly differ- 
ent values and the associated modes correspond to two 
different types of motions. 

In the limit of small Peclet numbers, the damping rate 
of the first type of mode is: 



These are exactly the weakly damped modes found in the 
context of the small-Peclet-number approximation. Note 
that, despite the high thermal diffusivity, temperature 
perturbations can be weakly damped, if they are associ- 
ated with vertical motions against the mean temperature 
gradient. 

The damping rate of the second type of mode is: 



Such modes are not found in the context of the small- 
Peclet-number equations. Note that this is not surprising 
since they correspond to solutions of the Boussinesq equa- 
tions which do not behave like Taylor series (see equations 
(||) and (||)) when the Peclet number goes to zero. These 
modes undergo a purely diffusive damping which can be 
made arbitrarily fast as the Peclet number vanishes. In- 
deed, whatever the values of k x and fc, , all these modes 



are reduced by an arbitrary large factor after a time pro- 
portional to P e /fc™ m . For this type of motions, the ver- 
tical advection term appearing in the linearized heat Eq. 
(||) is negligible. This shows that, in the limit of small 
Peclet number, temperature perturbations which are not 
produced by vertical advection are damped in a very short 
time. 

According to the above discussion, it is always possible 
to find a Peclet number such that, after an arbitrarily 
small time, the evolution of the infinitesimal perturbations 
is equally described by the Boussinesq equations or by the 
small-Peclet-number equations. 

We now consider another example of flow, yet in a lin- 
ear regime. It concerns the evolution of small disturbances 
in a stably stratified shear layer. This configuration dif- 
fers from the previous example by the presence of a mean 
horizontal flow sheared in the vertical direction. In this 
case, the validity of the linear version of the small-Peclet- 
number approximation has already been proved by Dudis' 
theoretical work (1974). This author considered specifi- 
cally an hyperbolic-tangent velocity profile in a stable at- 
mosphere characterized by a hyperbolic-tangent temper- 
ature profile and used a normal mode approach to study 
the stability of the flow. 

He first determined the neutral stability curve, i.e. 
the curve separating the stable and unstable regions in 
the parameter space, for decreasing values of the Peclet 
number. Then, he showed that for small Peclet numbers 
these neutral curves could be recovered using a linear 
version of the small-Peclet-number equations. The con- 
vergence of the Boussinesq equations towards the small- 
Peclet-number equations appears fairly rapid in this case 
since, already at P e — 0.2, the maximum difference be- 
tween the neutral curves is within 3 percent. We recently 
revisited the work of Dudis by considering a linear tem- 
perature profile instead of the tangent hyperbolic pro- 
file to characterize the stable stratification (Lignieres et 
al. 1999). We confirmed the validity of the small-Peclet- 



number equations to describe the neutral curves. In addi- 
tion, we verified its validity for other types of mode (un- 
stable modes symmetric to the shear layer mid-plane) as 
well as in the viscous case. Note that very rapidly damped 
modes corresponding to the second type of mode found in 
the previous discussion may also exist in this case. How- 
ever, they can not affect the stability of the shear layer 
since they are very strongly damped. 

The third example is a two-dimensional non-linear flow 
where a shear layer is forced at the top of a linearly 
stratified fluid. This flow has been studied numerically 
by Lignieres et al. ( 1998D for large Reynolds numbers 
(R e pa 2000 where R e is based on the layer thickness and 
the velocity difference across it). Figure la shows the typ- 
ical vorticity field resulting from the destabilization of the 
shear layer and the concentration of vorticity into vor- 
tices. The mean shear and the thermal stratification are 
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also represented in Figs, lb and lc (here, the means refer 
to horizontal averages). 

The other parameters being held fixed, we reduced the 
Peclet number from 1000 to 1 (cquivalently the Prandtl 
number P r has been decreased from 0.5 to 5 x 10~ 4 which 
already requires some computational effort). Figures Id 
and le present horizontal profiles of the vertical velocity 
and the temperature deviation for the two extreme values 
of the Peclet number. When this number is equal to 1000 
(Fig. Id), one recovers a classical property of the inflex- 
ional shear layer instability in stably stratified medium, 
namely that the phase lag between the vertical velocity 
and the temperature deviation is tt/2. By contrast, we 
observe that both fields are antiphased when the Peclet 
number is equal to unity (Fig. le). This striking differ- 
ence reveals a change in the predominant terms of the heat 
equation and we verified that this equation is now domi- 
nated by a balance between the vertical advection against 
the mean temperature gradient and the thermal diffusion. 
These first results are consistent with the convergence of 
the Boussinesq equation towards the small-Peclet-number 
approximation. A detailed comparison of the results of 
these simulations with those obtained by using the asymp- 
totic equations will be reported in a forthcoming paper. 

In this section, the validity of the small-Pcclct-numbcr 
approximation has been proved for two types of linear 
flows. For the case of infinitesimal perturbations in a lin- 
early stably stratified atmosphere, we noted that vertical 
length scales of perturbations have to be limited to fi- 
nite values to ensure uniform convergence. Moreover, the 
Boussinesq and the small-Peclet-number equations give 
the same solution after an arbitrarily small time, once ini- 
tial temperature perturbations not associated with verti- 
cal advection against the mean temperature gradient have 
been damped. The example of a non-linear flow we con- 
sidered is also consistent with the validity of the approxi- 
mation. 

4. Elementary properties of the 

small-Peclet-number approximation 

The elementary properties of the small-Peclet-number 
equations are analyzed in this section. From a practical 
point of view, the main interest of these equations is that 
their numerical integration does not require the compu- 
tation of the very rapid temporal variation of tempera- 
ture due to thermal diffusion. The Lagrangian derivative 
is indeed absent from the asymptotic heat Eq. (|l4|). This 
property is crucial for the investigation of small Peclet 
number regimes since numerical simulations are no longer 
limited by the huge separation between the dynamical and 
diffusive time scales. 

Another simplifying property appears when the Peclet- 
number equations (13), (|l4|), (p"5|), are written in terms of 



Using this rescaled temperature deviation, the small- 
Peclet-number equations become: 



du 



1 



— + u • Vu = -Vp + Rtpe z + — V z u, 



dt 

w = V 2 ^, 
V • u = 0, 
where 

R = R t Pe 



Re 



tr>t K 



(18) 

(19) 
(20) 

(21) 



Although these equations have been derived as a first or- 
der approximation of the Boussinesq equations in the limit 
of small Peclet number (see Sect. 2), they can also be in- 
terpreted as a zero order approximation of the Boussinesq 
equations in the limit of small Peclet number provided 
the non-dimensional number R — RiP e is assumed to 
remain finite. Starting from the dimensional Boussinesq 
equations, one only has to use P e AT, as a reference tem- 
perature instead of AT, and to assume Taylor-like expan- 
sions of the form (^) and . Then, provided R remains 
constant, the above equations arise at the zero order in 
P e . The derivation presented in Sect. 2 has been preferred 
because it does not require the assumption of an infinite 
Richardson number Ri. 

The main interest of the system @, (|^), (|o|), is that 
it only depends on two non-dimensional numbers, R and 
R e . This is an important simplification as compared to 
the original Boussinesq equations which are governed by 
three non-dimensional numbers, Ri, P e and R e . This sim- 
plification corresponds to the fact that the amplitude of 
the buoyancy force is no longer determined by two distinct 
processes, namely the vertical advection against the stable 
stratification which produces temperature deviations and 
the thermal diffusion which smoothes them out. It is now 
determined by a single physical process which combines 
the effects of both processes.. We shall show below that 
this process is purely dissipative and that this dissipation 
is anisotropic (not effective for horizontal motions) and 
faster for large scale motions. 

To do so, we write down the kinetic energy conserva- 
tion. Multiplying the momentum Eq. (|lj) by the velocity 
vector and integrating over the whole domain, we obtain: 



dt 



R w^idV 



1 

Re 



edV 



dS. (22) 



<V "-e Jv JS 

where, on the r.h.s of this equation, the first term is the 
work done by the buoyancy force, the second term repre- 
sents the viscous dissipation into heat and the third term 
is the kinetic energy flux on the surface bounding the do- 
main. Using Eq. (|19|), the work done by the buoyancy force 
can be divided into two terms to give: 

R I wifj dV = -R [ (V^) 2 dV + R [ VVV- • dS. (23) 
Jv Jv Js 
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Fig. 1. Effect of the Peclet number on a non- linear two-dimensional shear flow. The figures (a), (b), (c), show a 
typical vorticity field of the flow and the vertical profiles of the mean horizontal velocity and of the mean temperature, 
respectively. Horizontal variations of the vertical velocity and of the temperature deviation are displayed for two values 
of the Peclet number, P e = 1000 (d) and P e — 1 (e). Both quantities are antiphased in the low Peclet number case, 
this result pointing towards the validity of the small-Pcclct-numbcr approximation 



As temperature deviations vanish on the bounding plates, 
the second term also vanishes so that the kinetic energy 
conservation reduces to: 



dEi 



kin 



dt 



R I (Vip) 2 dV- 



1 

Rf, 



edV+ / F kin -dS.(24) 



This equation shows that the combined effect of the stable 
stratification and the thermal diffusivity is purely dissipa- 
tive. This simple result has to be compared with the case 
of the Boussinesq equations where the integrated work of 
the buoyancy could be positive or negative. As described 
in detail by Winters et al. (1995), it is then necessary 
to distinguish the amount of kinetic energy which is irre- 
versibly lost from the amount of kinetic energy which has 
been transformed into potential energy but can still return 
back to kinetic energy. Here, the situation is simpler since 
all the kinetic energy extracted by the buoyancy work is 
irreversibly lost. 



In order to specify the time scale of this dissipative 
process, we rewrite the above equations without the non- 
linear terms and for inviscid motions restricted to a ver- 
tical plane (e x ,e 2 ). The pressure term can first be elim- 
inated using the incompressibility condition (|20|). Then, 
the two momentum equations are combined to eliminate 
the horizontal velocity and the simplified heat equation 
allows to eliminate the rescaled temperature deviation. 
Finally, the evolution of the vertical velocity is governed 
by: 

dAAw 2 w 

-^T = R d^- (25) 

Considering isotropic motions of length scale I, the time 
scale of this process is l/(Rl 2 ), that is t\jt K in dimen- 
sional units. It appears that this dissipative process is 
faster at large scales than at small scales, which is just 
the opposite of what is observed in usual dissipative pro- 
cesses like thermal diffusion or viscous dissipation. Here, 
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however, the thermal diffusion does not act directly on 
the dynamics; it affects the temperature deviations which 
in turn modifies the buoyancy force amplitude. We have 
already seen that, in the limit of small Peclet numbers, 
rapid thermal exchanges lead instantaneously to a bal- 
ance between vertical advection against the mean strati- 
fication and thermal diffusion. This balance is described 



by equation (19) and it is straightforward to show that 



the resulting amplitude of the temperature deviations is 
stronger if the vertical velocity varies over a large length 
scale. The amplitude buoyancy force is therefore stronger 
for velocity fields varying over large length scales and this 
explains why the combined effect of the stable stratifica- 
tion and the thermal diffusion is faster at larger length 
scale. Note that, while classical dissipative processes are 
characterized by a Laplacian operator, the operator of the 
present dissipation is the inverse of a Laplacian. This can 
be seen by expressing ip as the inverse Laplacian of the 
vertical velocity and by reporting this expression in the 
momentum equation. 

Another interesting property is the anisotropy of this 
dissipative process. If one considers a velocity field of the 
form w cx exp[i(k x x + k z z)], where, as before, k z and k x 
represent its vertical and the horizontal scales, the char- 
acteristic time deduced from Eq. (Ea) is: 



kiy 



Rkl 



(26) 



which unsurprisingly corresponds to the inverse the damp- 
ing rate found in Sect. 2. The use of polar coordinates 
in the Fourier space is more appropriate to study the 
anisotropy of the process. With r 2 = fc 2 +fc 2 and tan(a) = 
kz/kx, the time scale becomes: 



1 



R \cos(a) 



(27) 



where a = tt/2 corresponds to horizontal motions and 
a = corresponds to vertical motions. 

We observe that the dissipation acts primarily on ver- 
tical motions while purely horizontal motions are not af- 
fected. This is not surprising since the buoyancy force 
only applies on the vertical component of the velocity. 
What is more interesting is that this anisotropy is stronger 
than in the context of the non-diffusive Boussinesq equa- 
tions. Indeed, for a given value of the wave vector mod- 
ulus, the time scale r increases faster towards horizon- 
tal motions (a — - > 7r/2) than the corresponding time 
scale of the buoyancy force in a non-diffusive atmosphere 
1/o\b = l/(cos(a)v r Ri). Considering motions strongly af- 
fected by the buoyancy force, we thus expect that these 
motions would be more predominantly horizontal in an at- 
mosphere dominated by thermal diffusion than in a non- 
diffusive atmosphere. 



5. Discussion 

In this paper, we derived a small-Peclet-number approx- 
imation, we discussed its validity for three flow examples 
and we analyzed its basic properties. In particular, we 
showed that the practical and theoretical difficulties char- 
acterizing the regime of very large thermal diffusivities 
and which had been mentioned in the introduction are 
considerably simplified in the context of the small-Peclet- 
number approximation. 

In what concerns applications for the dynamics of stel- 
lar radiative zones, it must be stressed that some type of 
motions can not be investigated using this approximation. 
First, there are no gravity waves in the context of the ap- 
proximation whereas these waves could play an important 
role in the radiative zone dynamics (Schatzman 1996). 
Second, thermal convective motions penetrating the ra- 
diative zone boundary have a high Peclet number so that 
the small-Peclet-approximation is not suitable to inves- 
tigate the overshooting layer at the boundary with the 
thermal convective zone. 

By contrast, there are various evidences that some mo- 
tions contributing to the radial transport of chemical ele- 
ments and angular momentum are characterized by very 
small Peclet numbers and could therefore be studied in the 
context of the small-Peclet-number approximation. The 
fact that the thermal structure is determined by the ra- 
diative heat flux only shows that the Peclet number char- 
acterizing eventual radial motions is necessarily smaller 
than unity. But other observational constraints, obtained 
by measuring the surface abundance of chemical elements, 
give much smaller Peclet numbers (see a recent review by 



Michaud & Zahn 1998). These Peclet numbers are defined 



as the ratio between the diffusion coefficient necessary to 
recover the observed surface abundance and the thermal 
diffusivity. In the absence of more sophisticated models, 
this diffusion coefficient is assumed to represent a verti- 
cal turbulent transport. For the sun, a Peclet number as 
small as 2 x 1CP 4 is obtained. This value may however be 
underestimated because the turbulence is most probably 
anisotropic. This as pect is taken into account in the model 
of Spiegel & Zahn ( 1992[ ) which describes the tachocline 
(the abrupt change in angular velocity at the top the so- 
lar radiative zone). The Peclet number which characterizes 
the horizontal turbulent motions and which is compatible 
with observed thickness of the tachocline remains much 
smaller than unity [P e sa 10 -2 ). 

These estimates suggest to use the small-Pcclct- 
numbcr approximation to investigate the property of small 
scale turbulent motions in stellar radiative zones. A first 
possible investigation could be the homogeneous turbu- 
lence in presence of a uniform mean shear and mean tem- 
perature gradient. With geophysical applications in mind, 
this configuration has already been extensively studied for 
large Peclet numbers (see Schumann 1996, for a review). 
A comparison with the small-Peclet-number case should 



8 



F. Lignieres: The small-Peclet-number approximation in stellar radiative zones 



be very instructive. Another important topic concerns the 
anisotropy between vertical and horizontal motions in an 
atmosphere dominated by thermal diffusion. Our linear 
study suggests that this anisotropy can be stronger than 
in a non-diffusive atmosphere. The ratio between vertical 
and horizontal turbulent viscosities could be affected and 
this can estimated through numerical simulations of the 
small-Peclet-number equations. 

Before concluding, it must be noted that the limit 
of lar ge diff usiviti es has already been considered (Spiegel 



1962| , Thual |1992[ ) in the context of the Rayleigh-Benard 



convection. However, an important physical property of 
thermal convection is lost in this limit. The thermal strat- 
ification is indeed assumed unchanged and this is not com- 
patible with the general observation that convective mo- 
tions transform the initial unstable stratification into an 
adiabatic stratification. There is no such an inconsistency 
for the case of stably stratified radiative zones we consid- 
ered here. 
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